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Abstract. 

We introduce and solve a special family of integrable interacting vertex models 
that generalizes the well known six-vertex model. In addition to the usual nearest- 
neighbor interactions among the vertices, there exist extra hard-core interactions 
among pair of vertices at larger distances. The associated row-to-row transfer matrices 
are diagonalized by using the recently introduced matrix product ansatz. Similarly as 
the relation of the six-vertex model with the XXZ quantum chain, the row-to-row 
transfer matrices of these knew models are also the generating functions of an infinite 
set of commuting conserved charges. Among these charges we identify the integrable 
generalization of the XXZ chain that contains hard-core exclusion interactions among 
the spins. These quantum chains already appeared in the literature. The present paper 
explains their integrability. 

Keywords: Integrable spin chains (vertex models), Solvable lattice models, Quantum 
integrability (Bethe Ansatz) 
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1. Introduction 

The six-vertex model was introduced by Pauling in order to explain the residual entropy 
of the ice at zero temperature. The model turns out to be of great interest for the physics 
and mathematics of many body interacting systems, due to its integrability. The row- 
to-row transfer matrix of the six- vertex model is the generating function for an infinite 
set of commuting non-trivial charges in involution pQ. The XXZ quantum chain being 
one of these charges. For a quantum system to be integrable its Hamiltonian should 
belong to an infinite set of commuting operators. The integrability of the XXZ quantum 
chain is then a consequence of the infinite number of commuting charges generated by 
the six-vertex model. For this reason the six-vertex model is considered as a paradigm 
of integrability in statistical mechanics |2j-[5j 

The connection of vertex models and quantum chains was also observed in a great 
variety of general vertex models and quantum chains thanks to the development of the 
quantum inverse scattering method [QISM] [6]-[8]. The QISM allows the construction 
of the vertex models and the associated quantum chains from the solutions of the Yang- 
Baxter equations (see [9j-[TT] for reviews). Basically we should expect that associated 
to any integrable quantum chain in one dimension should exists a vertex model whose 
row-to-row transfer matrix is the generating function for the quantum chain. 

On the other hand, along the recent years, it has been shown that it is possible to 
generalize several known integrable quantum chains by preserving their integrability. 
These generalizations are obtained through the introduction of suitable hard-core 
exclusion interactions. Examples of these generalizations were obtained for the XXZ 
quantum chain [T2J, [13] , spin-1 Fateev-Zamolodchikov and Izergin-Korepin models [14], 
SU(N) Sutherland and Perk-Schultz models [121 US] and the Bariev model [T6J, [T7] . 
Some of these extended quantum chains appear in the description of the stochastic 
dynamics of in the asymmetric exclusion problem with particles with extended sizes 
[T8| [15], [19] -[23]. The exact solution of these quantum chains were not obtained by the 
QISM but directly by the coordinate Bethe ansatz [24] or by the matrix product ansatz 
]25j-]28j. Consequently, the two dimensional vertex model generating these quantum 
chains are not known, and we cannot explain their integrability, up to now. We expect 
that similar to the quantum chains the associated vertex models should have, beyond 
the nearest neigbour interactions imposed by the lattice connectivity, also additional 
hard-core interactions. 

Recently an interacting five-vertex model in a square lattice was introduced and 
solved [29] . In this model there exist hard-core interactions along the diagonal of the 
square lattice. The exact solution was obtained only for the eigenspectra of the diagonal- 
to-diagonal transfer matrix, that distinctly from the row-to-row transfer matrix does not 
generate an infinite set of conserved charges. 

In this paper we are going to introduce a family of integrable interacting vertex 
models that generalize the six-vertex model and explains the integrability of the 
XXZ quantum chains with hard-core interactions. The row-to-row transfer matrices 
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associated to these vertex models generate an infinite set of conserved charges, that 
include the Hamiltonians of the hard-core interacting quantum chains. 

In the six- vertex model the interactions are between nearest neighbour vertices and 
are ruled by the geometrical connectivity of the allowed vertex configurations of the 
lattice. The interacting energy among two vertices is zero for allowed configurations 
or infinite otherwise. The family of interacting vertex models we introduce in this 
paper have additional interactions. Besides the usual nearest-neighbour interactions, 
imposed by the lattice connectivity, there exist also interactions among vertices at larger 
distances. 

The exact solution of transfer matrices associated to vertex models or quantum 
Hamiltonians are usually obtained through the Bethe ansatz [24] in its several 
formulations. The ansatz asserts that the amplitudes of the eigenfunctions of these 
operators are given by a sum of appropriate plane waves. Instead of making use of the 
Bethe ansatz, the solution we derive will be obtained through a matrix product ansatz 
(MPA). 

Under the general name of MPA several methods were introduced in the literature 
along the years. The first formulation was done for the description of the ground-state 
eigenfunction of some special nonintegrable quantum chains, the so called valence-bond 
solid models [30]-[33]. MPA becomes also a successfull tool for the exact calculation of 
the stationary probability distribution of some stochastic one dimensional systems [34J- 
[36J. An extension of this last MPA, called dynamical MPA was introduced in [371 [38] 
and extended in [39]. This last ansatz gives the time-dependent probability distribution 
for some integrable systems. 

The MPA we are going to use in this paper, in order to solve the new family 
of integrable vertex models, was introduced in [25]- [28]. This ansatz was applied with 
success in the evaluation of the eigenspectra of several integrable quantum Hamiltonians 
[25]-[27], transfer matrices [29], HQ] and the time-evolution operator of stochastic systems 
[28] . According to this ansatz, the amplitudes of the eigenfunctions are given in terms 
of a product of matrices where the matrices obey appropriated algebraic relations. In 
the case of the Bethe ansatz the spectral parameters and the amplitudes of the plane 
waves are fixed, apart from a normalization constant, by the eigenvalue equation of the 
Hamiltonian or transfer matrix. On the other hand, in the MPA the eigenvalue equation 
fixes the commutation relations of the matrices defining the ansatz. 

The layout of the paper is as follows. In section 2 we introduce the interacting 
vertex models we are going to solve. In section 3 we present the row-to-row transfer 
matrices associated to the vertex models. In section 4 we explain the MPA and give 
the exact solution of the vertex models. In section 5 we derive the associated quantum 
chains that commute with the row-to-row transfer matrices, solved in section 4. Finally 
in section 6 we conclude the paper with a general discussion. 
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2. The interacting four-vertex model 

The family of interacting vertex models we introduce are defined on a square lattice 
with M rows and L columns and toroidal boundary conditions. At each horizontal 
(vertical) lattice bond we attach an arrow pointing to the left or right (up or down) 
direction. These arrows configurations can be equivalently described by the vertex 
configurations of the lattice. A vertex configuration at a given site (center) is formed 
by the four arrows attached to its links. Similar to the asymmetric six-vertex model 
[47Tj |4~2] we impose that the allowed arrow configurations only contain vertices satisfying 
the ice rules: two of the arrows pointing inward and the other two pointing outward of 
its center. For the standard six-vertex model there are six possible configurations for 
the vertices. These configurations are showed in Fig. [T^l with their respective energies 
ei, . . . , Eq. The partition function is given by the sum of all possible vertex configuration 
with the Boltzmann weights given by the product of the fugacities e~^ 6i (i = 1, . . . , 6) 
of the vertices. In Fig. [lb, a more convenient notation is introduced, in which we only 
drawn, with the corresponding fugacities a , ai, b%, b 2 , C\ and c 2 , the arrows pointing to 
the left or down of the center defining the vertex. 




Figure 1. The six possible vertex configurations allowed by the ice rules, with their 
respective energies and fugacities. In (a) we draw all arrows and in (b) we give a 
representation where only the arrows pointing down and left are drawn. 



As a consequence of the ice rules and the periodicity of the lattice the number of 
down arrows n (or up arrows L — n) is conserved on each row of vertical bonds. 

It is important to notice that the ice rules imply the existence of interactions among 
the nearest-neighbor vertices. These interactions have a zero or infinite value. For 
example (see Fig. [lb), the vertex with fugacity oq has an infinite interaction energy if 
the vertex on its right side is one of the vertices (a 1; 61,02) or the vertex on its left side 
is one of the vertices (ai, &i,ci). The model with no extra interactions, besides those 
given by the connectivity of the vertices on the lattice, is the well known six-vertex 
model. The partition function of the model is given by the sum of all possible vertex 
configurations with the Boltzmann weights given by the product of the fugacities of the 
vertices. The six-vertex model is exactly integrable for arbitrary values of the fugacities, 
and is considered a prototype of an exact solvable model [2]- [4], [44], [42] . 
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We consider a special family of interacting four-vertex models. Besides the 
previously mentioned nearest-neighbor interactions (infinite or zero) imposed by the 
lattice connectivity, the models also contain interactions among pairs of vertices at larger 
distances. The allowed four vertex configurations, with their respective configurations 
are the vertices 1, 3, 5 and 6 shown in Fig. [TJ Contrary to the six- vertex model the 
vertex configurations with fugacities a\ and 6 2 are forbidden. Such interacting four- 
vertex models are labeled by two fixed positive integer S\ and s 2 that may take the 
values si, S2 = 1, 2, . . .. These parameters specify the additional interactions among the 
vertices that occur when there are vertices of types c\ and c 2 at distances equal to s\ 
or S2, in lattice units along the horizontal lines of the square lattice. A pair of vertices 
at distance (I = 1,2, ...), in units of lattice spacing, along this horizontal line interacts 
following the ruleajj 

(a) The interaction energy is zero if one of the vertices is a or b\. 

(b) If the pair of vertices is of type c% — c\ or c 2 — c 2 the interaction energy is infinity 
if I < si + s 2 . 

(c) If the pair of vertices is of type c\ — c 2 (c 2 — ci) the interaction energy is infinity if 
I < si (I < s 2 ) and zero if I > si (I > s 2 ). In the special case where I = s± (I — s 2 ) 
the interaction energy is finite and produces, besides the fugacity of the vertices, a 
contribution to the Boltzmann weight z\ (z 2 ) given by 

hb 2 , a ai . . 

Z\ = and z 2 = . (1) 

CiC 2 CiC 2 

The use of the missing fugacities a\ and b± to denote the interactions z\ and z 2 in (pQ) 
is convenient for the forthcoming analysis. In Fig. [2] we show some examples of allowed 
and not allowed configurations containing two arrows in the model with si = 2 ans 
82 = 1. 

In general, the contribution of a given pair of vertices is zero if the pair is not 
allowed (infinite interaction energy) or is given by the product of their fugacities (zero 
interaction energy). The exception to this rule happens when we have the pair c\ — c 2 
(c 2 — ci) at the distance s\ (s 2 ) along the horizontal line (see Fig. [2]). In this case, 
from (JTJ, the total contribution to the Boltzmann weight is given by C\C2Z\ = 6i& 2 
{c 1 c 2 z 2 = a ai). 

We can also extend our model to the cases where the interaction parameters s\, s 2 
may take the value zero. In these cases the interaction follows the previous rules (a)-(c) 
but instead of (CQ) we have the generalized form 

Z\ — Di ana Z2 — clq , I / ) 

CiC 2 " CiC 2 

where <5 Sj o is the usual Kronecker delta. In this generalization we may have the pair of 
vertices c\ — c 2 (c 2 — ci) at distance s\ — (s 2 = 0). On this case the total contribution 

I For simplicity we denote a given vertex by its fugacity, and a pair of vertices where the vertex fugacity 
ci is on the left (right) of C2 by c\ — C2 {c% — c±), respectively. 



Exactly solvable interacting vertex models 



(a) 



(b) 



2 2 ,4 

c 1 c 2 z 2 a b 



(c) 



(d) 



Figure 2. Configurations containing two arrows in a lattice of size L = 9, for the 
model with s± = 2 and s% = 1. The configurations (a) and (6) are allowed and 
their respective Boltzmann weights contributions are shown. Configurations (c) is not 
allowed due to the occurence of a pair c\ — c% at distance 1 < s\. configuration (d) is 
also not allowed because the two vertices c-i are at distance 3 < s\ + s-i- 



comming from the pair is now C\CiZ\ = 62 [c2c\Z2 = ax). The models where Sx = 
and S2 7^ (s 2 = and s\ 7^ 0) are equivalent to interacting five vertex models where 
only the vertex ax (62) of Fig. [T] is not allowed. The interactions on these models, 
according to rule (a), forbid arrows at distances smaller than s 2 + 1 (sx + 1). Moreover 
the pair of vertices c\ — c 2 (C2 — cx) interacts as in rule (c). The special case where 
Sx = S2 = recovers the standard six-vertex model, with interactions given only by the 
lattice connectivity and the ice rules. 

In summary, the general vertex models defined by the rules (a)-(c) and by (j2J), 
contains interacting four- vertex models (si 7^ 0, S2 7^ 0), interacting five- vertex models 
(sx = 0, S2 7^ or sx 7^ 0, S2 = 0) and the standard six- vertex model (sx = s 2 = 0). 

It is important to notice that the interactions, due to rule (b), forbid two vertical 
arrows at distances smaller than si + s 2 + 1 (sx, s 2 = 0,1,2,...). This can be interpreted 
as if the vertical arrows have an effective size sx + $2 + 1 (si, s 2 = 0, 1, 2, . . .), in units of 
lattice spacing. A vertical arrow on a given link has hard-core interactions that exclude 
the occupation of other vertical arrows at the link itself as well as the s± nearest links 
on its left and S2 nearest links on its right. In Fig. [3]we represent pictorically an allowed 
(a) and a not allowed (b) configuration of arrows for the interacting four-vertex model 
with Sx — 2 and s 2 = 1. In this example the arrows have an effective size sx + s 2 + 1 = 4. 
In the particular case sx = s 2 = the arrows have a unite size as it should be expected 
in the standard six-vertex model. 

Similarly as the connection of the six- vertex model with the XXZ quantum chain [4 J 
we are going to show in Sec. [3]that our interacting model, with general parameters sx and 
S2, is also related to quantum spin chains with hard-core exclusion interactions [121 E3]- 
The related models are the XXZ quantum chains with exclusion of up spins (a 2 -basis) 
at distances smaller than s\ + S2 + 1 [121 Ej and the time-evolution operator of the 
asymmetric exclusion problem of particles of size sx + s 2 + 1, in units of lattice spacing 

[El EH]. 
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Figure 3. Examples of an allowed (a) and not allowed (b) configuration of the 
interacting four- vertex model with s± — 2 and s% = 1. The vertical arrows have 
an effective hard-core size s = si + S2 + 1 = 4. 

3. The row-to-row transfer matrix of the interacting vertex models 

We denote a given configuration of n vertical arrows located at x\, . . . ,x n along the 
horizontal line of the lattice by \x\, . . . , x n ). The row-to-row transfer matrix T Sl)S2 , with 
elements 

(x\, . . . , x n \T sljS2 \x 1} . . . , x n ) = ag 1 6 1 3 c 1 5 c 2 6 2; 1 1 z 2 2 , (3) 

gives the Boltzmann weight contribution to the partition function due to vertical arrow 
configurations \xi, . . . , x n ) and \x[, . . . , x' n ) on two consecutive row^§|. In (J3J) are the 
numbers of vertices of type i (i — 1, 3, 5, 6) and n Zl (n Z2 ) is the number of pairs c\ — C2 
(c2 — ci) at distance of si (S2) units along the horizontal line (see Fig. [2]). The partition 
function of the interacting vertex model, due to the periodic boundary condition of the 
lattice, is given by 

Z = ^KJ- (4) 

As we discussed in the last section in the cases where s\ = 0, s 2 7^ and 
s i 7^ 0, S2 = the four-vertex models are equivalent to interacting five-vertex models. 
We can also verify directly that the corresponding transfer matrices To )S2 and T Sli o are 
related to those of models with si 7^ and s 2 7^ 0, namely, 

(b \ nsi / a \ nS2 

— J V Sl To iSl +s 2 — [j^J 'P S2 Ts 1 +s 2 ,0i (5) 

where V is the horizontal translation operator, i. e., V\x\, . . . , x n ) = \xi + l, . . . ,x n + l). 
This last relation enable us to relate transfer matrices of distinct models (si,S2) and 

§ In the particular case where we have no vertical arrows [n — 0) or in the case where s\ = S2 = 
(six- vertex model) and Xi = x[ (i = I,...,n) there exist two possibilities to connect the vertical 
arrows. In this case Eq. ^ has two terms and should be replaced by (x\, . . . , a; n |T SliS2 \x\, . . . , x n ) = 

niX'—n 1 l „L—n 
CLi0 1 + O2d 
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(s[, s' 2 ), provide si + s 2 = s[ + si 



2- 



r.x« = J = j 7> s ^t s , jS2 . (6) 

4. The Matrix Product Ansatz for the interacting vertex models 

The transfer matrix (j3J) with toroidal boundary conditions has a U(l) x Zf, symmetry 
due to the conservation of vertical arrows along the horizontal lines and the translation 
symmetry along the horizontal direction. Consequently the associated Hilbert space 
can be separated into block disjoint sectors labeled by the number n of vertical arrows 
(n = 0, 1, [L/(si + S2 + l)]) (we denote by [x] the integer part of x) and the momentum 
eigenvalues P (P = = 0, 1, ...,L — 1) of the translation operator V = e~ tP . We 

want to solve, in each of these sectors, the eigenvalue equation 

T sl>S2 \y niP ) = A„|* n ,p), (7) 

where A n and \^ n ,p) are the eigenvalues and eigenvectors of T SljS2 , respectively. Theses 
eigenvectors can be written in general as 

\^n,p) = A{x 1 ,...,X n )\x 1 ,...,X n ), (8) 

where A{x\, ...,£„) is the amplitude corresponding to the arrows configuration with n 
vertical arrows located at sites (xi, x n ), respectively. The symbol (*) in the sum 
means the restriction to the configurations where Xi+i — x% > si + s 2 (i — 1, . . . , n — 1), 
%i > 1) s\ — s 2 < x n — x\ < L — (si + s 2 ). Since \^ n ,p) is also an eigenvector with 
momentum P, then the amplitudes should satisfy 

■A-iXli Xn) _jp 



A{x\ + 1, ...,x n + 1) 
for x n < L, while for x„ = L we have 



e~ lP (9) 



A{x u ...,L) _ lV , 



A{\,x x + l,...,x n _i + 1) 6 ' ^ 

The exact solution of the transfer matrix ([3]) is obtained by an appropriate ansatz 
for the unknown amplitudes A{x\, x n ). As shown in the last section our model 
reduces to the standard asymmetric six-vertex model in the case where s\ = s 2 = 0. In 
this particular case the model is known to be exact solvable through the Bethe ansatz 
[2JH1JI12] an d also by the MPA [ID] introduced in [25]- [27]. In the present paper we are 
going to derive the exact solution of the interacting vertex models with s\, s 2 — 0, 1, . . . 
by using the PMA we propose in [27]. This ansatz asserts that to any amplitude in (jSJ) 
there exists a one-to-one correspondence with an ordered matrix product 

A(xx,...,x n ) E* 1 - 1 AE" 1 ' 1 A- ■ ■ £*"-*n-i-i AE L ~ Xn , (11) 

where the matrices A and E are associated with the sites where we have a vertical 
arrow or not, respectively. Actually the objects A and E are not necessarily matrices 
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but abstract operators with an associative producl[jj]. Their commutation relations will 
be fixed by imposing the validity of the eigenvalue equation of the transfer matrix 
(J7j). A well-defined eigenfunction is obtained, apart from a normalization factor, if all 
the amplitudes are uniquely related. Equivalently, in the subset of words (products of 
matrices) of the algebra containing n matrices A and (L — n) matrices E, should exist 
only a single independent word. The relation between any two words gives a c number 
that is equal to the ratio between the corresponding amplitudes in ffTTj) . 

Since the eigenfunctions \^ n ,p) (EJ) have a well-defined momentum P = ^-j 
(j = 0,...,L — 1), the relations (Q and (fit)]) imply the following constraints for the 
matrix products appearing in the ansatz (1111) 

E Xl ~ 1 AE X2 ~ Xl ~ 1 A ■ ■ ■ j7j Xn ~ x ^-i-~ 1 AE L ~ Xn = 

Q~ iP E Xl AE X2 ~ Xl ~ 1 A ■ ■ ■ jT 11 "" 1 "- 1 " 1 j^jjjL—xn—i (12) 

for x n < L, and for x n = L 

e -iP AE x x AE x*-xi-l A . . . E x n -x n -r-\ AE L-x n - X -\ _ ^3) 

In order to proceed, in the usual way, we are going to consider first the eigensectors 
of T ai)S2 with small values of n. 
The case n = 0. 

In this case the solution of the eigenvalue equation (J7|) is trivial since we do not 
have vertical arrows between two successive rows. There are only two possible horizontal 
arrangements: either all bonds have a horizontal arrow or all of them are empty. In this 
case the vertices are all of type 1 or type 3 (see figure 1) and consequently the eigenvalue 
is given by 

A = 4 + &f. (14) 
The case n = 1. 

We have in this case just one vertical arrow between two rows. The transfer matrix 
links a vertical solid line at position x (x = 1, . . . , L) above a row to a vertical line at 
any position y (L > x > y + s\ or y — S2 > x > 1) under this row. The elements of the 
transfer matrix T sl S2 (y,x) in this sector with momentum P are given by ([3]). They are 
given by the product of the Boltzmann weights of the vertex appearing on the row. The 
vertex configuration at the sites x and y will be of types C2 and c\. If the position of the 
line x is greater or equal (less or equal) than y and all the others vertices will be of types 
b\ (a ) and ao (&i) depending on whether the vertices are between the positions x and y, 
or not, respectively. However in the special cases where x = y + Sj or x = y — s 2 we have 
also to include the fugacities Z\ (z 2 ) due to the long range interaction (j2j) between these 
vertices. Consequently the eigenvalue equation (CO) for the transfer matrix (jHJ) associated 

| In the original formulation [25} [26] of the ansatz an unnecessary additional matrix [43] and a trace 
operation was taken in the right hand side of (jTTJ) . The MPA in (fTT]) is formulated as in the general 
formulation given in |27] . 
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with the components of \ip n .p) fill I) with n = 1 and momentum P, give us the algebraic 
relations for the matrix products 

x—si—l 

k x E x - x AE L - x = o^ x ^ l b x - y - x c x c 2 E y ' x AE L - y 

y=l 

L 

+ E al- x - x b\- y+x - x c x c 2 E y - x AE L ^ 

y=x+S2 + l 

+ a^-^b^hE^^AE 1 --^ 31 + b^-^a^a^+^AE 1 ^- 32 , (15) 

where we inserted (J5J) in the last two terms in the right hand side. 

The equation (fl5j) is simplified by expressing all the matrix products in terms of a 
single one. This is done by imposing that the matrix A depend on a spectral parameter 
k. Without loss of generality, the matrix A can be written in terms of the matrix E and 
a spectral parameter dependent matrix Ajjl 

A = E^AtE 1 ' 32 , (16) 
with Ak satisfying the following algebraic relation with the matrix E 

EA k = e tk A k E. (17) 
Inserting (TIB"]) in ({15]) and using (fTTI) we can factorize (|T5|) : 

{Ax - X jr 1 a^~ x+y ~ 1 bt~ v - 1 c 1 c 2 e ik ^ - £ af^ 1 ^-^ 1 

y=l y=x+s 2 + l 

x Cl c 2 e ik{y - x) - a^^b^he-^ 31 

- b^-^a^a^jE^^AkE 1 -^ 1 - 32 = 0. (18) 

In order to produce a nonzero norm state we should impose E x ~ 3l ~ l AkE L ~ x+1 ~ 3 ' 2 ^ 0, 
for x = 1, L. As a consequence we must have, by evaluating the sums in ([181) . 

A x = a L L{k) (h.ye- lk3l +b^M(k) (^e*' 



where 



r / 7 \ a &2 + (cic 2 - 6ife 2 )e ifc , M a ai - cic 2 - aihie lk 

L y k > = 2 r~ ik ' M fc = 1 — mk • 20 

— a 0ie lk a bi — bfe %K 

In order to satisfy (fl5j) . the eigenvalue Ai in (fl9|) should independ on the vertical line 
position x. Thus the last term in the right hand side of (fl9|) must vanish. The only way 
to cancel this term, for non zero Boltzmann weights, is obtained by imposing e lkL = 1, 
that fixes the spectral parameter k 

k=^J (j = 0,l,...,L-l). (21) 

The most general relation A = ^2™ =1 A a Akj E@ could be used. However (fTI>)) is more convenient 
since otherwise the structure constants S(kj, fa) that will appear in (|32p will depend on the parameters 
Si and S2- 
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The eigenvalue (1191) is then given by 



Ai = a^L(k) f ^ j e" ifcsi + b^M(k) (j-) ' e ikS2 , (22) 



where the values of the spectral parameter k are given by (J2TJ). 

Finally, by inserting ([IS]) and using (ITT)) in the equations (IT2l) and (IT3)) we verify 
that coincides with the momentum of the eigenstate: 

k = P=— 3 (j = 0,l,...,L-l). (23) 
The case n = 2. 

In this sector there are two vertical arrows in the row. We have in general two 
types of relations for the amplitudes (ITT)) . Relations where at least one of the vertical 
arrows (1/1,1/2) are at distance si or S2 from (x\,x 2 ), and those where y\ and y 2 interlace 
with x\ and x 2 . (1 < y\ < X\ — s\\ x\ + s 2 < y 2 < x 2 — Si) and (x\ + s 2 < y\ < x 2 — s%; 
x 2 + s 2 < y 2 < L) . Then, the eigenvalue equation (JTj) imply 

X\—S\ X2—S1 * 

A 2 E X1 - 1 AE X *- X1 - 1 AE L ~ X * = E aS~ (x2 ~ Vl+1) C2 

2/1=1 V2=Xl+S2 

x ^(xx, yi )/(y 2 , xi)»(x 2> y 2 )EK- 1 AE"-K- 1 AE L -y* 

X2—S1 L * 

+ E E b L 1 -^ +1) c 1 f(y 1 ,x 1 ) 

x <?(x 2 , yi )f(y 2 , x 2 )E^- x AE^- x AE L -y\ (24) 

where the symbol * in the sums means that terms with \y 2 — yi\ < Si + s 2 are excluded 
and 

, / c 2 al~ x ~ X if y > x + s 2 
[ a 2 ai/ci if y = sc + s 2 

/ s / b\~ y ' x c x if x > y + s x 
I o 1 1 6 2 /c 2 if x = y + si 



The relation (124)) connects configurations where the arrangements of vertical arrows 
(x\,x 2 ) above the row do not have the same distance of the vertical arrows (y\,y 2 ) 
bellow the row. In other words, the distance of the incoming arrows y 2 — y\ are in 
general different of the outcoming distance x 2 — x\. To solve the relation ( 124")) we need 
now to use a generalization of the algebraic relation (ITBI) for the case of two arrows. 
The generalization of ( IT6j) is done by writing the matrix A in terms of two new spectral 
parameter matrices as: 

A = j2E-^A k] E 1 ~ s \ (27) 
j=i 

with the commutation relations 

EA kj = e^A k .E (j = l,2), (28) 
where the spectral parameters k\ and k 2 are up to now unknown complex numbers. 
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Inserting (1271) in (1241) and using (I2"gj) and (|2"U|) we obtain, after similar manipulation 
as we did in the case n — 1, the following constraint 



2 

E 

x M{kj)M{k{) }e ik i Sl e- iklS2 e ik i Xl e ihX2 A k A kl 



E« 

i,z=i 



L(ki)M(kj) 



aib 2 



V 
o 



S1-S2 



3,1=1 

x e i ( fc j+ fc i)( s l- s 2) e «(fcj+fc;)x 2 ^ fc ^ 



L{ki)M{ki) 



aib 2 
a bi 



22— 



V 

a . 



;,i(fej+fci)(si-S2) 

~ s 2 / 7 \ 1 

\ao/ 



+ E4 

'&iY 



r 2 2 
c l c 2 



6f(a - 6ie-* fc i)(a - &ie~ ifc 



x 



V 



-E 



3,1=1 



r 2 2 



i(oo - &ie ifc O(«o - he- ik ') 



Xi 



\a J \a 



+ CJo C l C 2 a l 
3,1=1 



giLki gifcj(s2+si) p -ik t S2 p ikjX\ 



6f (a — &ie _ 



q^ikj si ^ikixi 

b\ (do — b\e~ tk i 



6i 
a 



.l"2 



+ E fl C l C 2 & 2 
3,1=1 

a o / V a o y 



QiLki e -ikiS2 e ikjX2 pi/ysip— ifti(s2+si)pjfcz:C2 



X 



ao&i(ao — feie lfci ) ao&i(ao — ^ie" 
^fc^fc, = 0, 



(29) 



where 1 < Xi < X2 < L with |x2 — xi| > S\ + S2 + 1. This can only be satisfied if each 
sum is identically zero. Moreover since A2 should be independent of x\ or x 2 a possible 
solution of (J2"9~|) is obtained by imposing 

2si 



A 2 = a r 



h \ / \ 2s2 

^ e-^+^LikMh) + &f ( - ) e^ 1 ^* 2 
ao 7 V bi J 



\aoJ 
x M(k 1 )M(k 2 ). 



(30) 



The algebraic relation between the matrices A kl and Afc 2 are obtained by imposing that 
both the second and third terms in ff29l) are zero independently, i. e., 

A kj A kl — S(kj,ki)A kl A kj , (A kj f = (j ^ Z = 1,2), (31) 
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where 

L{kj)M(ki) - ^ 

S(k 3 M = :- m)M{kj) (32) 

with and M(fc) given by fl20l) . Note that the structure constants S(kj,ki) is 

independent of si and S2 due to the choice fl27|) (see the footnote related to Eq. ffTB"]) ). 
Finally, the vanishing of the last four terms in ( 129]) will give us relations that fix the 
spectral parameters values k\ and /c2- These equations are obtained by exploring the 
algebraic relations fl3TT) 



Cik \ s i+ s 2 
J^j Sfo,*,) (I, J = 1,2 and Z^j). (33) 

The eigenvalues are obtained by inserting the solutions (ki,k 2 ) of these last equations 
in (J30l) . The momentum P is obtained by using (1271) and ( J28l) in (fl2l) and ( TTBl : 
P = ki + k 2 . The amplitudes of the corresponding eigenvectors are obtained from 
the algebraic relations (|27D. fl28|) and (1311). 
The case of general n. 

The general case follows straightforwardly from the n = 2 case. The previous 
calculation can be extended for arbitrary values of the number n of vertical arrows. The 
eigenvalue equation for the transfer matrix connects the amplitudes A(xi, . . . ,x n ) and 
A(yi, ...,y n ) where: 1 <yi< x 1 -s 1 ,xi+s 2 <yi< x 2 -s x , ■ ■ . ,x n _i+s 2 <y n < x n -s 1 
and x\ + S2 < yi < x 2 — si, £2 + ^2 < Z/2 < ^3 — si, . . . , x n — s 2 < y n < L. To solve the 
eigenvalue equation we need to extend the definition (|27j) and the commutation relations 
(J2HJ) for general n, i. e., 



A = Y,E- Sl A k] E 1 - s * (34) 

with 

EA k . = e ik *A kj E (j = l,...,n). (35) 

The parameters kj (j = 1, . . . , n) are in general unknown complex numbers that will be 
fixed by the eigenvalue equation ([71). Inserting ( 1341) in (J7J) and using the commutation 
relations (|3~5"I) we obtain, similarly as done in the case n = 2, the algebraic relations 
among the matrices A k . (j = 1, . . . , n): 

A k] A kl = S{kj, h)A kl A kj (A k] ) 2 = (j^l = l,...,n), (36) 

where S{kj,k{) is given by (|32l) and the spectral parameters fcj (j = l,...,n) are fixed 
by the equation 

e iL kj = _ JJ ^_ j k) {j = 1; ^ n) _ (37) 

The acceptable set {/cj} of spectral parameters defining the eigenvectors \^ n ,p) is given 
by the solutions of (I3TI) where A;; 7^ fcj (j, / = l,...,n). Since (^4^) = 0, solutions of 
with coinciding roots give us null states. 
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No new algebraic relations appear for the matrices {^4fc-} besides ( 13 5 p and (13T)j) . 
This not easy to verify. The difficulty here is similar to that in showing the absence 
of many-particle scattering in the standard coordinate Bethe ansatz for vertex models 
[2]- [5]. The product associativity of the algebra (135]) and (|36|) follows from the property 
S(kj, ki)S(ki, kj) = 1. From the the eigenvalue equation (JTj) we obtain the eigenvalues 

C, \ nsi n \ns2 n 

-1 J] e-^ Sl L(%) + b L A f) J] e^Mik,), (38) 
a oJ j=1 \bxJ j=l 

where L(k) and M(k) are given by (120]) and the spectral parameters {kj} are the 
solutions of (13"7I) . Finally, the momentum P follows from f[T2"j) . ffTSl . f[3~4"|) and (1331) : 

n 

P = £*i- (39) 
i=i 

Before closing this section let us give a possible representation for the matrices E 
and A of the ansatz (ITT]) [27]. For a given solution {fei, k n } of the spectral parameter 
equations fl37|) . in the sector with n vertical arrows, the matrices E and {A^, ...,ylfc n } 
have the following 2 n x 2 n dimensional representation: 

E = ®( 1 n L), (40) 




•v (8)1 " ; I ®( ;; ) ; j, (41) 

where S(kj, ki) are given by (1321) and A is obtained from (1341) . It is important to notice 
that the matrices E and Au have the same functional form obtained for the formulation 
of the MPA for the Hamiltonian of the XXZ spin chain [27]. It happens because the 
matrices E and satisfy the same algebraic relations (I3~4"|) - (I3T)|) in both models. The 
only difference are the values of the spectral parameters {kx, k n } that for the present 
model are fixed by equation (1371) . 

It is important to notice that all the interacting four- vertex (si 7^ 0, S2 7^ 0), five- 
vertex (si = 0, S2 7^ and si 7^ 0, S2 = 0) and six-vertex models (si = S2 = 0) have the 
same functional form for the operators A^. (j = 1, . . . , n). The only difference being the 
structure constants S(kj,ki), since the spectral parameters ki,...,k n are solutions of 
distinct equations. As a test of our calculations we see that (|37l) . for the case si = S2 = 
coincides with the spectral parameter equations obtained by the well known Bethe ansatz 
calculation of the asymmetric six- vertex model [HI H2]. Using the representation (|40l) 
and (HT1) we also obtain on this case the same amplitudes for the eigenfunctions as those 
derived by the Bethe ansatz. 

We can also solve the interacting models introduced in this paper by the use of the 
same coordinate Bethe ansatz used in the six- vertex model [2]- [5]. In this case instead 
of obtaining the structure constants fl32l) . with no dependence on s\ or S2, we obtain 
a S-matrix that depends on these parameters. Although the final spectral parameter 
equations (1371) and the eigenfunctions for both ansatz are the same, we believe the 
derivation through our MPA is more elegant. 
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The phase diagram of the interacting vertex models introduced in this paper can 
be derived from that of the six-vertex model. This follows from the spectral-parameter 
equations ( 1371) . that can be written as 

n 

e i[L-n( S1 + S2 )]k je iP(s 1+ s 2 ) = _ JJ s ^ y = 1; ^ n ^ ^ 

1=1 

where we have used fl39l) . Consequently the eigenvalues belonging to the eigensector 
with n vertical arrows and momentum P of the transfer matrix T Sl>S2 with parameters 
si and S2 are related to those of the standard six- vertex model (s\ = S2 = 0). The 
related six-vertex model is defined on a cylinder with perimeter L' = L — n(si + S2) and 
with a momentum-dependent seam spanning its lengthEJ 



5. The interacting vertex models and generalized XXZ quantum chains 

It is well known the connections among the XXZ quantum chain and the standard 
asymmetric six- vertex model [I] . The row-to- row transfer matrix of the six- vertex model 
is the generating function of an infinite set of commuting charges. The XXZ quantum 
Hamiltonian being one of these charges. Since the six-vertex model is a particular case 
(si = S2 = 0) of our general model it is natural to expect that we can also generate 
generalized exact solvable XXZ quantum chains for the general case where s\, S2 > 0. 

The first step towards the derivation of those generalized quantum chains is the 
identification of the row-to- row transfer matrix T SljS2 , for arbitrary values of s\ and S2, 
clS cl generating function of commuting charges. 

It is important to notice that the structure constants S(kj, h) given in ( 1321) do not 
depend on the particular values of s\ and s%. It only depends on the parameters a , a\, 
b%, b 2 , ci and cq,. Consequently it is convenient to use the same parametrization used in 
[H] for the asymmetric six-vertex model 



a = e"^ 1 


_ ch+v smh(j + en) 
sinh7 




ai = e- fe 


_ c _ /l _ i ,sinh(7 + e/i) 
sinh7 




bi = e"^ 3 


_ e _ h+l/ sinh/i 
sinh7 


(43) 


b 2 = e"^ 4 


_ ^.^sinh/i 
sinli7 




ci = e"^ 5 


= 1 




c 2 = e"^ 6 


= 1, 




where e = ±1 and 7, /i 


, h and v are free parameters. 


The choice C\ — — 1 does 



not restrict the model since, due to the periodic boundary condition, these parameters 

+ The phase e lP( - Sl+s ^ in the left hand side of jUJ can be obtained by considering a vertex model 
with twisted boundaries. The twisted is obtained by introducing a seam with distinct vertex fugacities 
along the vertical direction. 
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always appear in pairs. The symmetric interacting models are obtained by setting 
h = v = 0. 

Inserting in ( 1321) the parameters in ( T43j) we obtain 



_ 1 + e 4h e i ^ +k ^ - ee 2h cos je ik < 

^ J ' l ' ~ \ _|_ e 4/i e i(fcj+fc() _ eg 2h CO S7e ifc -'' ' 

that shows the independence of the structure constants of the algebra defining the MPA 
on the parameters fi and v. 

The results of previous sections imply that, for a given solution {hi, . . . , k n } of the 
spectral parameter equations (1571) . the amplitudes of the corresponding eigenfunction 
of T S1>S2 depends only on the algebraic relations (1M|) - (!36|) . Consequently all the 
eigenfunctions of T SljS2 , for fixed values of s\ and S2, are the same for arbitrary values 
of fi and v, i. e., 

[T sl , S2 (n,v),Ts 1 , S2 (n',v')} =0. (45) 
Moreover, since from (j^J) T SlyS2 and T s > tS > differ by the multiplication of a diagonal 
operator, we may write the general relation 

[T SliS2 (/i, v), T s > i>s > 2 (/j,' , u')] =0, si + s 2 = s[ + s' 2 - (46) 

The conserved charges are derived by expanding T SuS2 ([m,v) in terms of \i. Since 
we are interested only in the derivation of the quantum chain associated to T S1)S2 (/^, u), 
in the following we are going to consider only the two leading terms in the /z-expansion. 

The matrix elements of T S1)S2 (//, v) in the eigensector with n arrows are given bjQ 

({y}\T SuS2 (^u)\{x}) = P s ({y})P s ({x})a L - {x "- yi+1) c 2 

x g(xi, yi)f(V2, xi)g(x 2 , 2/2) • • • f(y n , x n -i)g(x n , y n ), (47) 

if x\ > yi + si, and 

({y}\T sljS2 M\{x}) = P s ({y})P s ({x})bi-^ +1) c 2 

f(yi,x 1 )g(x2,yi)f(y 2 ,x 2 ) ■ • -g(,x n ,y n -i)f(y n ,x n ), (48) 

if Xi < y\ — s 2 . In these last equations s = Si + s 2 and P s ({x}) projects out configurations 
not satisfying the hard-core exclusion constraint: 

n-l 

P s ({x}) = 9(x n — Xi — [L — (s + 1)]) J] 6(x i+ i - Xi -(s + 1)), (49) 

i=l 

where 

*>-{;£ :^ 

is the standard step function. In fj47j) and fHHj) the functions /($/, 2) and g(x, y) are 
defined in ( 1251) and ( 1261) and have the leading behaviour in the /x-expansion 

f(y,x) = ^ + u)W- x - Sl -2) {1 + y _ x _ Si _ 2) + 2S y ,_ XiS+1 ] eyUC oth 7 } 
x %' - x - (s + l))^y- si -i + O^ 2 ), (51) 

* The exceptional cases considered in the footnote of page 7 should be considered separately, with no 
changes in the final results. 
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-.u—h 



g(x,y) 



sinh7 



x—y',0 



+ 



.y—h 



sinh7 



x—y',1 



+ e 



h—v i 



>x-y'-\/ 



x <W-*i-i + 0(/i Sl+2 ) 



(52) 



Using ([BT]) and ([32]) in fllTJ we obtain 

({y}|T sl)S2 |{x}) = 



sinh7 



x 1 + 



sinh 7 

2 E^Si-^+i) ecosh 
i=i / 



0(// 2 ))p s ({*}), 



7 



(53) 



where 



"li+l— Xj,S + l 



(54) 



^i+i-xi.s+l for i = l,...,n-l 
v 5 Xn _ Xli L-( s +i) for i = n 

We can rewrite (153]) in terms of spin-| Pauli matrices by identifying an arrow at site % 
as of = +1, and we obtain 

Y^L crf+1 
/ \ Si > 



- Zjj=l 2 



sinh 7 



sinh 7 / 
H sliS2 +0{^) 



V 



(55) 



where 



E 

,i=i 



+ r 



-/i m ^(i + <)}p si 
i=i j 



+ S2 5 



e 2/ \ e_ 



e _2h h 



Si 



A = ecosh 7, e = ±1. 



(56) 



In this last equation periodic boundary conditions is imposed, of = (o x ±ia y )/2 and 
P Sl + S2 is the diagonal operator ({x}\P Sl+S2 \{x}) = P Sl+S2 ({x}) defined in (SHI), that 
projects out configurations where two up spins are distances smaller than si + s 2 + 1. 
The same expressions (155"]) and (T5B1) are also obtained by inserting ( I5T1) and ( l52l) in ( )48l) . 

We see from ( 156]) that associated to the family of solvable interacting vertex 
models, with interacting parameters Si and s 2 , we have a generalized asymmetric XXZ 
quantum chain. This quantum chain have hard-core exclusion interactions that forbid 
the occupation of up spins at distances smaller than si+s 2 -|-l. The important parameter 
controlling the exclusion effect is given by the combination s = s± + s 2 . The effect of 
the parameter si in (156]) is the same as an external magnetic field h m — 

The Hamiltonian (156]) for the case where s\ = s 2 = gives the standard asymmetric 
XXZ chain. The particular case where A = e + + e_ gives the time-evolution operator 
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of the asymmetric exclusion problem where the particles have hard-core size si + S2 + 1, 
in units of lattice spacing [TU [15] . 

The integr ability of ( |56|) . for general values of s\ and s 2 , is already known in the 
literature [TU [151 EH] • What is unknown is why the inclusion of the hard-core exclusion 
constraints in the standard XXZ quantum chain model (si = S2 = 0), producing the 
Hamiltonian fl56|) . does not destroy the integrability of the original model. This paper 
gives the explanation for this fact. The quantum chains fl56|) . for a given value of 
and s 2 , belongs to the infinite set of of commuting charges generated by the row-to- row 
transfer matrix T Sl+S2 of our interacting vertex models. 

6. Conclusions and generalizations 

We have introduced in this paper a special family of exact solvable interacting vertex 
models. These models are generalizations of the six-vertex model. Besides the usual 
nearest-neighbor interactions imposed by the lattice connectivity, the models also 
contains hard-core interactions along the horizontal lines of the lattice. The range of the 
additional interactions depends on two fixed integer parameters Si,s 2 = 0,1,.... This 
new family of models contain interaction four- vertex models (si 7^ 0, S2 7^ 0), interacting 
five-vertex models [s\ = 0, s 2 7^ 0, or s\ 7^ 0, s 2 = 0) and the standard six-vertex model 
(s% = S2 = 0). These vertex models can also be interpreted as if the vertical arrows 
entering in the vertex configurations have an effective hard-core size s = si + s 2 + 1. 

For a given value of s\ and s 2 the model can be parametrized by four parameters 
(7, z/, /i and h) (see ( H3|) ). The exact solution of the eigenspectra of the row-to-row 
transfer matrix shows the independence of the eigenfunctions on the parameters v and /x, 
implying the commutativity of an infinite family of transfer matrices (see ( 1451) and ()46l) ). 
The expansion of the transfer matrix in terms of the free parameter /x give an infinite 
set of commuting conserved charges. A member of this set is a XXZ quantum chain 
with hard-core interactions that exclude two up spins (<r 2 -basis) at distances smaller 
than s\ + s 2 + 1 [121 [13] . This quantum chain also describes the asymmetric diffusion 
of particles with hard-core size s\ + s 2 + 1, in units of lattice spacing [TBI [15], [19]. The 
integrability of this quantum chain is already known [121 EE] ■ The present paper give an 
explanation for this integrability, namely, the existence of an infinite set of commuting 
charges generated by the row-to-row transfer matrix of the introduced interacting vertex 
models. 

The phase diagram of the introduced interacting vertex models, as discussed in Sec. 
4, can be obtained from the known results of the six-vertex model. The critical phases 
of the models are governed by a Coulomb gas type of conformal field theories with 
central charge c — 1. Since the spectral parameter equations are the same as those of 
the extended XXZ quantum chain [13] the compactification ratio, that fixes the critical 
exponents, depend on the density of vertical arrows n and on the parameters s\ and s 2 . 

The exact solution we obtained for the interacting vertex models introduced in this 
paper were done by using the matrix product ansatz introduced in [25j-[28j. The most 
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elegant solution of the six-vertex model is obtained by using the i?-matrix approach in 
the quantum inverse scattering method. Can we solve the introduced interacting vertex 
models by generalizing the R matrix of the six-vertex model? Due to the appearance 
of nonlocal operators in the row-to-row transfer matrix we were not able to find such 
generalization. That is an interesting question to be answered in the future. 

Similar to the hard-core generalization of the XXZ quantum chain, has been 
also generalized several known integrable quantum chains through the introduction 
of suitable hard-core exclusion interactions. As examples we have hard-core 
generalizations of the spin-1 Fateev-Zamolodchikov and Izergin-Korepin models [14|, 
S q U(N) Sutherland and Perk-Schultz models [19], [2^ and the Bariev model (TTJ [16]. 
These quantum chains with no hard-core interactions, as the standard XXZ chain, are 
known to be generated by the row-to-row transfer matrix of known vertex models. We 
believe that the same ideas presented in this paper for the generalization of the six- 
vertex model can also be extended to these vertex models. We expect that the vertex 
models generating the quantum chains with only next neighbour interactions are special 
cases of a larger family of integrable interacting vertex models. 
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